Self-avoiding walks on scale-free networks 
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Several kinds of walks on complex networks are currently used to analyze search and navigation in 
different systems. Many analytical and computational results are known for random walks on such 
networks. Self-avoiding walks (SAWs) are expected to be more suitable than unrestricted random 
walks to explore various kinds of real-life networks. Here we study long-range properties of random 
SAWs on scale-free networks, characterized by a degree distribution P(k) ~ k~~' . In the limit of 
large networks (system size N — > oo), the average number s n of SAWs starting from a generic site 
increases as //", with fi = (k 2 ) / (k) — 1. For finite N, s n is reduced due to the presence of loops in 
the network, which causes the emergence of attrition of the paths. For kinetic growth walks, the 
average maximum length, (L) , increases as a power of the system size: (L) ~ N a , with an exponent 
a increasing as the parameter 7 is raised. We discuss the dependence of a on the minimum allowed 
degree in the network. A similar power-law dependence is found for the mean self-intersection length 
of non-reversal random walks. Simulation results support our approximate analytical calculations. 

PACS numbers: 89.75.Fb, 87.23. Ge, 05.40.Fb, 89.75.Da 
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I. INTRODUCTION 

Many natural and artificial systems have a network 
structure, where nodes represent typical system units 
and edges represent interactions between connected pairs 
of units. Thus, complex networks are currently used to 
model several kinds of real-life systems (social, biologi- 
cal, technological, economic), and to study different pro- 
cesses taking place on them [1-3]. In recent years, new 
models of complex networks have been designed to ex- 
plain empirical data in several fields. This is the case 
of the so-called small-world [4] and scale-free networks 
[5], which incorporate various aspects of real systems. 
These complex networks provide us with the underlying 
topological structure to analyze processes such as spread 
of infections [6,7], signal propagation [4,8], and random 
spreading of information [9,10]. They have been also em- 
ployed to study statistical physical problems as percola- 
tion [6,11] and cooperative phenomena [12,13]. 

In a scale-free (SF) network the degree distribution 
P(k), where k is the number of links connected to a 
node, has a power-law decay P{k) ~ k~~< . This kind 
of networks have been found in social systems [14], for 
protein interactions [15], in the internet [16], and in the 
world-wide web [17]. In both natural and artificial net- 
works, the exponent 7 controlling the degree distribution 
is usually in the range 2 < 7 < 3 [3,18]. The origin 
of such power-law degree distributions was addressed by 
Barabasi and Albert [5], who argued that two ingredi- 
ents are sufficient to explain the scale-free character of 
many real-life networks, namely: growth and preferen- 
tial attachment. They found that the combination of 
both criteria yields non-equilibrium SF networks with an 
exponent 7 = 3. One can also study equilibrium SF 
networks, defined as statistical ensembles of random net- 
works with a given degree distribution P(k) ~ k^ 1 [3], 



for which one can analyze several properties as a func- 
tion of the exponent 7. SF networks display the so-called 
small- world effect, and they have been found to be ultra- 
small, in the sense that the mean distance between sites 
increases with the network size N slower than log N [19]. 

Social networks form the substrate where dynamical 
processes such as disease propagation and information 
spreading take place. These networks have the property 
of being searchable, i.e. people (nodes in a network) 
can direct messages through their network of acquain- 
tances to reach a distant specific target in only a few 
steps [20-22]. It is clear that the structure of such net- 
works will play an important role in these dynamical pro- 
cesses, which are usually studied by means of stochastic 
dynamics and random walks. Several characteristics of 
random walks on complex networks have been analyzed 
in connection with diffusion and exploration processes 
[23-25]. In this context, it is known that some processes, 
such as navigation and exploratory behavior are neither 
purely random nor totally deterministic [26] , and can be 
described by walks on graphs [27,28]. 

Self-avoiding walks (SAWs) can be more effective than 
unrestricted random walks in exploring a network, since 
they cannot return to sites already visited. This prop- 
erty has been used by Adamic et al. [28] to define local 
search strategies in scale-free networks. However, the 
self-avoiding property causes attrition of the paths, in 
the sense that a large fraction of paths generated in a 
stochastic manner have to be abandoned because they 
are overlapping. This can be a serious limitation to ex- 
plore networks with pure SAWs. 

SAWs have been traditionally used to model structural 
and dynamical properties of macromolecules [29,30]. 
They are also useful to characterize complex crystal 
structures [31] and to study critical phenomena in lattice 
models [32]. Universal constants for SAWs on regular 
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lattices have been discussed by Privman et al. [33]. In 
our context of complex networks, the asymptotic prop- 
erties of SAWs have been studied recently in small- world 
networks [34]. 

Here we study long-range properties of SAWs on equi- 
librium scale-free networks, and discuss the 'attrition 
problem'. The number of surviving walks to a given 
length n is obtained by an approximate analytical proce- 
dure, and the results are compared with those obtained 
from numerical simulations. In particular, we find that 
the number of surviving walks after n steps scales as 
a power of the system size N. We note that the term 
'length' is employed throughout this paper to indicate 
the (dimensionless) number of steps of a walk, as usually 
in the literature on networks [3]. 

The paper is organized as follows. In Sec. II we give 
some definitions and concepts related to SAWs, along 
with details on our computational method. In Sec. Ill wc 
calculate the number and end-to-end separation of SAWs 
in (uncorrelated) scale-free networks. In Sec. IV we an- 
alyze the length at which non-reversal random walks 
intersect themselves in these networks (self-intersection 
length), and in Sec. V we calculate the average attrition 
length of kinetic growth SAWs, at which they cannot 
continue without violating the self-avoidance condition. 
The paper closes with some conclusions in Sec. VI. 



II. BASIC DEFINITIONS AND METHOD 



A self-avoiding walk (SAW) is defined as a walk along 
the bonds of a given network which can never intersect 
itself. The walk is restricted to moving to a nearest- 
neighbor site during each step, and the self-avoiding con- 
dition constrains the walk to occupy only sites which have 
not been previously visited in the same walk. 

The simplest procedure to obtain SAWs consists just 
in generating ordinary random walks and stop when they 
arrive at a node already visited. A problem with this 
sampling algorithm in regular lattices is the exponen- 
tially rapid attrition for long walks, since the probability 
of an n-step walk being self-avoiding behaves for large 
n as e~ A ™, where A is the so-called attrition constant 
[35]. Due to this limitation, more sophisticated schemes 
based on Monte Carlo sampling, have been employed to 
generate SAWs with the correct weight, and to obtain 
ensemble averages of several quantities [35] . This has al- 
lowed, for example, to model the equilibrium statistics of 
linear polymers in dilute solutions. In general, for net- 
works including nodes with different degrees (contrary to 
usual regular lattices), sampling by using simple random 
walks introduces a bias in the weight of different SAWs. 




FIG. 1. Schematic diagram showing a non-reversal random 
walk of length n = 5 on a realization of a random graph. 
Open and black circles represent unvisited and visited nodes, 
respectively. The starting node is indicated by a larger circle. 
The non-reversal condition allows in principle for the next 
(sixth) step three possible nodes (denoted A, B, and C). The 
self-avoiding condition excludes node C for the sixth step. For 
a non-reversal SAW one chooses among nodes A, B, and C. 
If C is selected, then the walk stops. For a kinetic growth 
walk, one chooses A or B. 

One can also consider kinetically grown SAWs, which 
can be more adequate to analyze dynamic processes. 
Such walks are well-suited to study, for example, search 
or navigation processes on networks, where they are as- 
sumed to grow step by step in a temporal sequence. In 
the following we will consider two kinds of growing walks. 
The first kind will be 'non-reversal' self-avoiding walks 
[35]. In these walks one randomly chooses the next step 
from among the neighboring nodes, excluding the previ- 
ous one. If it happens that one chooses an already visited 
node, then the walk stops (see Fig. 1). These walks will 
allow us to study the 'self-intersection length' (see Sec. 
IV). The second kind of walks considered here are kinetic 
growth walks [36], in which one randomly chooses the 
next step among the neighboring unvisited sites and stops 
growing when none are available. These walks were stud- 
ied to describe the irreversible growth of linear polymers 
[36], and will allow us to consider the 'attrition length' 
for a walk on a given network (see Sec. V). Note that 
kinetic growth walks are less sensitive to attrition than 
non-reversal SAWs, since in the former the walker always 
escapes whenever a way exists. 

We consider SF networks with degree distribution 
P(k) <~ fc~ 7 . They are characterized, apart from the 
exponent 7 and the system size TV, by the minimum de- 
gree fco, which affects markedly some characteristics of 
SAWs in these networks (see below). We assume that 
P(k) = for k < ko. Our networks are uncorrelated, in 
the sense that degrees of nearest neighbors are statisti- 
cally independent. This means that the joint probability 
P(k, k') fulfills the relation [3] 

P(k,k') = ^P(k)P(k'). (1) 

For the numerical simulations we have generated net- 
works with several values of 7, k 0} and N. To generate 
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a network, once denned the number of nodes Nk with 
degree fc, we ascribe a degree to each node according to 
the set {Nk}, and then connect at random ends of links 
(giving a total of L = ^2 k kNk/2 connections), with the 
conditions: (i) no two nodes can have more than one 
bond connecting them, and (ii) no node can be connected 
by a link to itself. We have checked that networks gener- 
ated in this way are uncorrelated, i.e. they fulfill Eq. (1). 
All networks considered here contain a single component, 
i.e. any node in a network can be reached from any other 
node by traveling through a finite number of links. For 
each set of parameters (7, ko, N), we considered different 
network realizations, and for a given network we selected 
at random the starting nodes for the SAWs. For each 
considered parameter set, the total number of generated 
SAWs amounted to about 5 x 10 5 . 

For regular lattices, the number s n of different SAWs 
starting from a generic site has an asymptotic depen- 
dence for large n [33]: s n <~ n r_1 /i n , where T is a critical 
exponent which depends on the lattice dimension, and fj, 
is the so-called 'connective constant' or effective coordi- 
nation number of the considered lattice [37] . In general, 
for a lattice with connectivity ko, one has fj, < ko — 1. 
This parameter \x can be obtained as the limit 



u = lim 



ri— >oo s 



71-1 



(2) 



The connective constant depends upon the particular 
topology of each lattice, and has been determined very 
accurately for two- and three-dimensional lattices [35]. 

For poissonian and scale-free networks the number of 
SAWs of length n depends on the considered starting 
node of the network. In the sequel we will call s n the 
average number of SAWs of length n, i.e. the mean value 
obtained (for each n) by averaging over the network sites 
and over different network realizations (for given 7, fco, 
and N). For Erdos-Renyi random networks with poisso- 
nian distribution of degrees, one has s^f — (k) n [34], and 
therefore the connective constant is fi = (k). In connec- 
tion with this, we note that for a Bethe lattice (or Cayley 
tree) with connectivity fco, the number of SAWs is given 
by Sn L = &o(fco — l) n ~\ an d one has ubl = ko — 1. 



III. GENERAL CHARACTERISTICS OF SAWS 

A. Number of walks 

We calculate first the average number, r„ , of different 
n-step unrestricted walks starting from a node chosen at 
random. One trivially has n = (fc). To calculate r n for 
n > 1, one needs the degree distribution for nodes at 
which one arrives following a random edge. Thus, given 
a generic node and a link starting on it, we call Q(k) the 
degree distribution for the other end of the link. The 
probability of reaching a node with connectivity k is pro- 
portional to fc; therefore 



(3) 



where (fc) in the denominator is a normalization factor. 
Then, the average number of two-step random walks is 
given by r 2 = ri{k)q, where the subscript Q indicates 
that the average value is taken with the probability distri- 
bution Q(k). We find r 2 = (fc 2 ) [average values without 
subscripts are taken with the degree distribution P(k)]. 
For n > 2 we have r n — r„_i(fc)g, and then 



(fc) 



(fc 2 ] 
(fc) 



(4) 



We have checked that this expression for r n and those 
given below for SAWs coincide with those derived by 
using a generating function for the degree distribution 
[28,38]. 

We now calculate the average number, s n , of differ- 
ent self-avoiding walks of n steps starting from a node 
taken at random, in uncorrelated networks. We will first 
consider the case n/N — > (thermodynamic limit). For 
n = 1, one has si = (fc). For n > 1 we take into account 
that each (n — l)-step walk arriving at a node with de- 
gree fc, gives rise to fc — 1 n-step walks. Thus, we have 
s n = s n -i(k — 1)q, which yields 



(fc 2 ) 
(fc) 



Then, the connective constant // 
by 



(5) 



for TV — > 00 is given 



f-^oc 



(fc 2 ) 
(fc) 



- 1 . 



(6) 



This is consistent with the fact that uncorrelated net- 
works arc locally tree-like, and the ratio (fc 2 )/(fc) is the 
average degree of a randomly-chosen end node of a ran- 
domly chosen link [3]. For SF networks with 7 < 3, fi 
diverges as N — > 00, as a consequence of the divergence of 
(fc 2 ). For 7 > 3, we can approximate the average values 
in Eq. (6) by replacing sums by integrals, and find 



fco 



7-2 
7-3 



1 . 



(7) 



Note that the ratio s„/s„_i does not depend on n for 
system size N — > 00. This is equivalent to assume that 
nodes in different steps of a non-reversal random walk 
are different. This means, in other words, that the prob- 
ability of finding loops with n' < n in a n-step walk is 
negligible. For finite networks, however, there will appear 
loops of any size [39] , introducing corrections to the num- 
ber of SAWs, and s n will be lower than given by Eq. 5. 
These corrections will be of order n/N for n/N <C 1. The 
effects of this reduction in the number of random SAWs 
in scale-free networks will be considered in Sections IV 
and V. 
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As indicated above, the number of SAWs on regular 
lattices scales for large nass„~ n r_1 /i n , where T is a 
critical exponent which depends on the lattice dimension 
D, and one has T = 1 for D > 4 [33,35]. For the SF 
networks studied here we find s n <~ (U^, indicating that 
T = 1, the same exponent as for regular lattices in many 
dimensions. 



B. End-to-end separation 

For walks on regular lattices, one usually considers 
an end-to-end Euclidean distance. Our SF networks, 
however, lack a metric and a true distance is not de- 
fined. Thus, we will consider the end-to-end separation 
for SAWs on SF networks as a function of the walk length 
n, the separation between two nodes being the number 
of links along the shortest path connecting them. In Fig. 
2 we present the average separation (d n ) from the n'th 
node in a SAW to the starting one (n = 0) for SF net- 
works with 7 = 3, k = 3, and several sizes N. This 
average separation (d n ) increases first linearly with n 
and finally saturates to a finite value that depends on 
the system size. This is logical if one takes into account 
that these networks are locally tree-like, and for the first 
steps the minimum separation between nodes n and is 
d n = n. As n increases, there appear shorter ways con- 
necting nodes n and 0, and finally (d n ) becomes inde- 
pendent of n. This saturation for relatively small values 
of n is consistent with the small average separation be- 
tween pairs of nodes (diameter) in this kind of networks. 
In fact, it is known that the diameter increases slowly as 
In N/ In In ./V for 7 = 3, and even more slowly as In In A 
for 7 < 3 [19]. The limit of (d n ) for large n in our SAWs 
is lower than the diameter of the considered networks, 
because nodes with larger k have a higher probability 
of being visited than those with smaller k, as indicated 
above. Since the average separation of a node to all other 
nodes in a network decreases for increasing degree of the 
selected node, the average separation (d n ) for large n in 
SAWs is lower than the diameter of the network. For ex- 
ample, for two of the system sizes represented in Fig. 2 
(N = 2.1 x 10 3 and 7.7 x 10 4 ) we have for large n, (d n ) = 
4.60 and 5.88 vs diameters of 4.77 and 6.07, respectively. 

For regular lattices, the mean squared end-to-end dis- 
tance of SAWs scales for large n as n 2v 1 v being a 
dimension-dependent critical exponent. The upper criti- 
cal dimension for these walks is D = 4, [40] which means 
that above this dimension one has v — i, as for brow- 
nian motion (Markovian random walks). On the other 
side, for the SF networks studied here, the mean squared 
end-to-end separation scales as n 2 in the thermodynamic 
limit, i.e., with an exponent v = 1. This exponent co- 
incides with that corresponding to D = 1, reflecting the 
fact that loops become irrelevant in the considered net- 
works as N — > oo (networks become tree-like). This be- 
havior is, however, not captured by SAWs on finite SF 



networks, for which {d n ) converges to a constant, and 
therefore gives a null exponent for sufficiently large n. 
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FIG. 2. Average separation between the n'th node and 
the starting one, for kinetically-grown SAWs on scale-free 
networks with 7 = 3 and fco = 3. Symbols corre- 
spond to different system sizes. From bottom to top: 
N = 2.1 x 10 3 ,9.6 x 10 3 ,7.7 x 10 4 and 6.2 x 10 s . Dotted 
lines are guides to the eye. 



IV. SELF-INTERSECTION LENGTH 

As indicated above, finite-size effects on SAWs on any 
(finite) network will be appreciable as soon as the walks 
are long enough, as a consequence of the presence of loops 
in the network. Thus, SAWs are a suitable tool to probe 
the large-scale topological structure of complex networks. 
In particular, the probability of a walk intersecting itself 
will depend on the system size, as well as on the topol- 
ogy of the network under consideration. To study this 
probability, we consider here non-reversal self-avoiding 
walks [35], that stop when they try to visit a node al- 
ready visited in the same walk. The number of steps 
of a given walk before intersecting itself will be called 
'self-intersection length' of the walk, and will be denoted 
I. 

To obtain the mean self-intersection length (l) of these 
walks, we will calculate the probability that a walk stops 
at step n (<ti N). Let us consider for the moment nodes 
with a given degree k. The average number of nodes with 
degree k visited after n steps is 

V k = n Q(k) , (8) 

and the average number of those yet unvisited is 



U k = N k - V k = NP(k) - nQ(k) . 



(9) 



Then, the probability of reaching in step n an unvisited 
node is u k oc kUk, and that of finding one already vis- 
ited is Vk oc (k — 2)Vk- This is due to the fact that a 
visited node has k — 2 possible links to reach it, as two of 
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its connections are not available because they were em- 
ployed earlier: one for an incoming step and one for an 
outgoing step. Therefore, the probability p n of finding in 
step n a visited node with any degree is 



EfeK + u k) ' 



(10) 



Inserting into Eq. (10) expressions (8) and (9) for V4 
and Uk, and keeping terms linear in n/N one has 



Pn 



(k)N 



]T(fc-2)Q(fc), 



and finally 



where 



n 

Pn^W- 



W = 



(fc 2 )-2(fc) 
(fc> 2 



(11) 



(12) 



(13) 



Note that for the networks considered here w > 0. In 
fact, (k 2 ) — 2(k) > is the condition to have a giant 
component in a network [41]. 
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FIG. 3. Fraction of non-reversal SAWs that survive after 
n steps, without intersecting themselves. Results are plotted 
for SF networks with 7 = 3, fco = 3, and several system sizes 
N. From left to right: N = 3.3 x 10 3 , 2.6 x 10 4 , 7.7 x 10 4 , and 
2.1 x 10 5 . Solid and dashed lines indicate results of numerical 
simulations and analytical calculations, respectively. 

To calculate the probability distribution for the self- 
intersection length I, we consider Mq random walks start- 
ing from nodes taken at random. We call M\ (n) the num- 
ber of non-reversal SAWs that remain after n steps (i.e., 
those which did not find any node visited earlier). Thus, 
M\(n) — M\(n + 1) = p„Mi(n), and considering n as a 
continuous variable x, we have a differential equation for 
Mx(x): 



1 dM-L x 
M^~chT ~ ~ W N 

which yields for integer n: 



Mi (n) = M exp 



w n 
' "2 7V 



(14) 



(15) 



In Fig. 3 we present results for the fraction of surviv- 
ing walks Mi(n)/M for SF networks with 7 = 3. We 
compare the curves derived from Eq. (15) (dashed lines) 
with those obtained from numerical simulations for net- 
works with fco = 3 (solid lines). Both sets of results agree 
with one another in the sensitivity region of our numeri- 
cal procedure (down to <~ 10~ 6 ). For larger n (and lower 
Mi(n)/M ) one expects terms of higher order in n/N to 
become relevant, and Eq. (15) to be less reliable. 
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FIG. 4. Mean self-intersection length (l) as a function of 
system size, for SF networks with fco = 3 and several values 
of the exponent 7. From top to bottom: 7 = 5, 3, 2.5, and 
2. Symbols are results of numerical simulations and dashed 
lines were obtained from Eq. (17). Error bars of simulation 
results are less than the symbol size. 

The average self-intersection length of these walks can 
be obtained as (I), with the probability distribution 



fiffl I ( wl 2 
*(0 = «,-exp(--- 



(16) 



which gives the probability of returning to a visited site 
in step I. Treating I as a continuous variable, we replace 
sums by integrals and find 



(I) 



nN 
~2w 



(17) 



For large N and 7 > 2, (k) converges to a finite value, 
and w ~ (k 2 ); therefore, the mean self- intersection length 
scales as (I) - {N/{k 2 )) 1 / 2 . For 7 > 3, (k 2 ) does not 



5 



diverge for large N, and then (/) ~ VAT. Mean self- 
intersection lengths derived from Eq. (17) are shown 
in Fig. 4, along with those found in numerical simula- 
tions for several exponents 7. Both methods give results 
agreeing with one another within the error bars of the 
numerical simulations. In general, we find (I) <~ N 13 ', 
with an exponent j3 that decreases from 0.5 to 0.25 as 7 
decreases from 3 to 2. From the distribution for Z given 
by Eq. (16) one finds a mean-square deviation for the 
self-intersection length of the walks: of = CN/w, with a 
constant C = 2 — tt/2. This means that a;/(Z) w 0.52. 

We note that the probability distribution R(l) (and 
therefore the average value (/)) is independent of the min- 
imum degree fco- It depends, apart from the system size, 
on the exponent 7 of the degree distribution through the 
mean values (fc) and (fc 2 ). In Fig. 5 we present the mean 
self-intersection length for SF networks with 7 = 2.5 and 
several values of ko, as derived from numerical simula- 
tions. Different Ze values give a unique A^-dependence 
of (I). Thus, the dependence of (Z) on fc , which should 
appear in Eq. (17) through the dependence of w on ko, 
is negligible for our purposes. However, the minimum 
degree ko affects strongly other properties of SAWs, such 
as the attrition length studied in the following section. 
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FIG. 5. Mean self-intersection length (l) as a function of 
system size for SF networks with 7 = 2.5 and several values 
of the minimum degree fco: squares, ko = 3; circles, fco = 5; 
triangles, fco = 7. 



V. ATTRITION LENGTH 

In this section we consider random SAWs, that travel 
on the network until they arrive at a node (called here- 
after 'blocking node'), where they cannot continue be- 
cause all adjacent nodes have been already visited and 
are not available for the walk. These are kinetic growth 
walks, as defined by Majid et al. [36]. The number of 
steps of a given walk until being blocked will be called 
attrition length of the walk, and will be denoted L. 



We now calculate the average attrition length of kinetic 
growth walks, and obtain its asymptotic dependence for 
large system size N. To find this average length we will 
derive a probability distribution for L, in a way similar 
to that employed above for the self-intersection length. 
With this purpose we note that a blocking node for a 
kinetic growth walk is characterized by the fact that all 
its links except one (employed for an incoming step) con- 
nect it with nodes previously visited. Then, for step n of 
a walk and for a given degree k, the average number N' k 
of blocking nodes present in the network is given by the 
binomial formula: 



N' k = kN k p k n -\l- Pn ) 



kN k p k n ~ 



(18) 



where p n (<C 1 for n -C N) is the average fraction of links 
joining a generic node with nodes visited earlier, as given 
in Eq. (12). This means that N' k ~ fc 1_7 p£ _1 , and then 
one has N' k <C N' k for k > ko. Thus, the average number 
of links connecting the n'th node in a walk with nodes 
already visited coincides, within our approximation, with 
N' kQ . (Note that there is one such link available for each 
of those N' ko nodes.) 

The probability of finding a blocking node in step n is 
given by the ratio q n = J^k N' k /N en d, where N en d = (k)N 
is the total number of ends of links in the network. Then, 
to order n/N, we have q n w N' kg /(k)N. This approxima- 
tion relies on the fact that p n = wn/N -C 1, which may 
be unfulfilled when the minimum degree ko is large and 
the average attrition length can be on the order of the 
system size N (see below). 

To derive the probability distribution for the attrition 
length L, we consider Mo kinetic growth walks starting 
from nodes taken at random. We call M 2 (n) the number 
of walks that survive after n steps. Then, the number of 
walks finishing at step n (for which n is a blocking one) 
is M 2 (n) — M 2 (n+1) = <7„M 2 (n). Considering again n as 
a continuous variable x, one has a differential equation 
for M 2 (x): 



1 dM 2 



= -Y 



(5) 



fco-l 



M 2 dx 

with the network-dependent constant 

^0 fc ...fcp-1 

N (k) 



(19) 



(20) 



and w given in Eq. (13). Then, for integer n we have 

fco" 



M 2 (n) = Mo exp 



(21) 



which gives the number of walks that remain after n 
steps, i.e. M 2 (n)/Mo is the probability of surviving to 
length n. Here x is a number (dimensionless length) 
given by x k Q a = k Q N k °- 1 /Y. 

Therefore, the probability distribution Z(L) for the at- 
trition length of these walks is 
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-L\ ko 



(22) 



This distribution is strongly dependent on the minimum 
degree fco, since nodes with this degree are in fact con- 
trolling the maximum length of kinetic growth walks in 
scale-free networks. The distribution Z(L) is displayed 
in Fig. 6 for 7 = 3, fco = 3, and different system sizes. 
Dashed lines were obtained from Eq. (22), and solid lines 
were derived from numerical simulations. Both sets of 
results follow the same trend, but the numerical results 
seem to be larger than the analytical ones for large L. 
This difference is larger than the noise of the simulation 
results, and shows the validity limit of our approximation 
for large L. 



degree (fc). For 7 < 3, we have (L) <~ N a , with an expo- 
nent a that decreases from 1 — l/fc to (1 — l/fco)/2 as 7 
decreases from 3 to 2. 

In Fig. 7 we show the average attrition length (L) 
as a function of the system size N for 7 = 3. Symbols 
correspond to results of numerical simulations for several 
values of the minimum degree fco, and dashed lines were 
obtained by using Eq. (23). For the largest fco, (L) de- 
rived from simulations increases with N slightly faster 
than the analytical result. This difference is not strange 
if one observes that for fc = 9, (L) is on the order of TV 
(in fact, for N — 10 5 , N/(L) w 3), and our assumption 
that n <C N for all steps of SAWs is not true. However, 
even in this case Eq. (23) gives a rather good approxi- 
mation for the average length (L) (see Fig. 7). 



N 




12000 



FIG. 6. Distribution probability Z(L) for the attrition 
length L of kinetic growth walks on scale-free networks with 
7 = 3, fco = 3, and different system sizes. From left to right: 
N = 3.3 x 10 3 ,2.6 x 10 4 ,7.7 x 10 4 , and 2.1 x 10 s . Solid and 
dashed lines indicate results of numerical simulations and an- 
alytical calculations, respectively. 
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FIG. 7. Mean attrition length (L) as a function of system 
size for SF networks with 7 = 3 and several values of fco- 
From top to bottom: fco = 9, 5, and 3. Symbols are results 
of numerical simulations, and dashed lines were obtained by 
using Eq. (23) . Error bars of simulation results are less than 
the symbol size. 



From the distribution Z(L) we obtain an average at- 
trition length 

r being Euler's gamma function. Thus, the dependence 
of (L) on N for large systems is controlled by xq. To 
obtain the asymptotic dependence of Xo, we note that 
Y in Eq. (20) scales for 7 > 2 as Y <~ w ko ~ 1 , because 
Nk /N converges to a constant for large N. In addi- 
tion, w - (fc 2 ), and therefore x h Q ° - (N "/ '(fc 2 ))^" 1 . For 
7 > 3, (fc 2 ) converges to a finite value as N — > 00, and 
the average attrition length increases as (L) ~ _/V 1_1//fc °. 
For 7 = 3, (fc 2 ) - In AT, and (L) - (N/ lnTV) 1 - 1 /^ This 
means that for a given system size, the average number of 
nodes visited in kinetic growth walks rises with increas- 
ing fc , as a consequence of the increase in the average 



In order to define strategies to search in this kind of 
networks, it is clear that nodes with low degree limit 
the effectiveness of the process. Thus, actual strategies 
based on SAWs have to include additional conditions to 
improve their efficiency. In this line, Adamic et al. [28] 
have proposed an algorithm based on SAWs that prefer 
high-degree nodes to low-degree ones. In any case, the 
long-range properties of pure SAWs give us direct insight 
into the structure of SF networks, further than the local 
neighborhood of a node, where the structure of links is 
tree-like. On a larger scale, one always finds loops sooner 
or later in finite networks, which is in fact probed by 
SAWs. In particular, the average self-intersection length 
(I) given in Eq. (17) is a measure of the typical size 
of loops in equilibrium SF networks. The presence of 
loops in a network is responsible for attrition of the walks. 
Then, the mean attrition length (L) given in Eq. (23) is 
a measure of the long-range 'openness' of a network (the 
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longer (L), the less loops contains a network). In this 
sense, our SF networks become more 'open' as their size 
N increases, and eventually are loop-free (or tree-like) in 
the thermodynamic limit, where both (I) and (L) diverge 
to infinity. 

As a result, we find that the efficiency of SAWs to 
explore scale-free networks increases for increasing expo- 
nent 7. This is a consequence of the fact that for a given 
system size N, the fraction of nodes with high degree in- 
creases for decreasing 7. High-degree nodes are visited 
more probably than low-degree ones, and once visited the 
former are more effective to block a SAW in later steps 
(they have more connections), thus reducing the mean 
self- intersection and attrition lengths. 



VI. CONCLUSIONS 



Self-avoiding walks provide us with an adequate tool 
to study the long-range characteristics of SF networks. 
In particular, they allow us to study the quality of a net- 
work to be explored without returning to sites already 
visited. For large networks, the number of SAWs in- 
creases as s n /s n -\ w (k 2 )/(k) — 1, provided that n <C N. 
For a given n, s n decreases with decreasing system size, 
as a consequence of the presence of loops in the net- 
works. These finite-size effects affect strongly the max- 
imum length of kinetic growth walks on scale-free net- 
works. 

We have calculated self-intersection and attrition 
lengths by using an approximate probabilistic method, 
which yields results in good agreement with those de- 
rived from numerical simulations. Both, the average self- 
intersection length and attrition length scale as a power 
of the system size N. For the mean self-intersection 
length of non-reversal SAWs we have (I) ~ with j3 
depending on the exponent 7 of the degree distribution. 
In particular, for 7 > 3 one has j3 — 0.5, and decreases as 
7 is lowered. The length of kinetic growth walks in scale- 
free networks is limited by attrition of the paths, and the 
mean attrition length follows a dependence (L) <~ N a , 
with a depending on 7 and the minimum degree k . For 
7 > 3, one has a = 1 — 1/fco- This dependence of the 
exponent a on fco is remarkable, reflecting the fact that 
the length of SAWs is limited by attrition at sites with 
the minimum degree fc . 
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